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Abstract 

We consider how, for quasi-degenerate neutrinos with tri-bi-maximal mixing at a high-energy scale, the 
mixing angles are affected by radiative running from high to low-energy scales in a supersymmetric theory. 
The limits on the high-energy scale that follow from consistency with the observed mixing are determined. 
We construct a model in which a non-Abelian discrete family symmetry leads both to a quasi-degenerate 
neutrino mass spectrum and to near tri-bi-maximal mixing. 

1 Introduction 

Neutrino oscillation data is at present consistent jl], ^ with just three light neutrinos with near tri-bi-maximal 
(TBM) mixing between flavours |3|, ||, [|^| [r|. However the nature of the mass spectrum is still not established, 
being consistent with cither a normal or an inverted hierarchy. Moreover, although the magnitude of the mass 
squared difference between neutrinos is reasonably well determined, the absolute scale of mass is not, being 
consistent with both a strongly hierarchical spectrum or a quasi-degenerate (QD) spectrum. 

Radiative running is especially important for QD neutrinos, as the effects on mixing angles are larger for 
QD neutrinos than in the hierarchical case. This was stressed in B, ^ where the mixing favoured at the time, 
bi-maximal mixing, was studied in-depth. More recent studies of mixing angles running include |l^, O, [l^, |l3) 
(and references therein). Here we discuss radiative corrections to TBM mixing, assuming that it arises through 
new physics, such as a family symmetry, at a high-energy scale. We determine how high, in a supersymmetric 
extension of the Standard Model, the initial energy scale can be while maintaining near TBM mixing at the 
low-energy scales relevant to oscillation experiments. The main difference from existing work is that emphasis 
is placed on the energy scales rather than on the resulting low-energy angles. Specifically, we set the angles 
to their TBM values at high-energy scales, run the angles to low-energy and iterate the process to find the 
highest-energy scale that still keeps the low-energy angles within current experimental bounds. The process 
is then repeated for different points of the parameter space, and the results are presented as a contour plot in 
the m Vi — tan/3 plane (i = 1 for normal and i = 3 for inverted hierarchy). 

The underlying question raised by the observed near TBM mixing is the origin of the pattern and the reason 
it is so different from quark mixing. Models based on family symmetries, particularly discrete non-Abelian 



family symmetries, have been constructed to explain this pattern, e.g. [[14, 15 1. In these models the difference 
between the quark and lepton sector follows naturally from the see-saw mechanism together with a strongly 
hierarchical right-handed neutrino Majorana mass spectrum. However these models only apply to the case of 
an hierarchical neutrino mass spectrum. Here we discuss how a discrete non-Abelian family symmetry can 
also give rise to near TBM mixing for the case of a QD spectrum. 



2 Radiative corrections to TBM mixing 

Family-symmetry models are typically constructed at some high scale, Mf, at which the model specifies 
relationships among parameters. To compare the predictions to low-energy data, radiative effects should be 



* ivo@cftp.ist.utl.pt 

tg.ross@physics.ox.ac.uk 

tserna@physics.ox.ac.uk 



1 



Normal Hierarchy Inverted Hierarchy 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 



m vl (eV) m v3 (eV) 

Figure 1: Shows contours of Log 10 (Afp) where Mp is the highest-energy family-symmetry breaking scale at 
which we can set TBM and have the neutrino mixing within 4cr of the low-energy observed values. The white 
region in the lower left of the contour plots are the regions where Mp can be greater than 10 16 GeV. 

considered through the use of the renormalization group equations. When there is a strong hierarchy, it is 
often the case that these running effects do not change the mixing angles by much jn], |ll], |l^, |l3| . In the case 
of QD neutrinos, however, the mixing angles can change a lot with the energy scale, to the point of erasing any 
special structure arranged by a family symmetry. For model-building purposes it is very important to know 
the highest-energy scale at which we can start with TBM mixing and still be consistent with mixing-angle 
data after running the angles down to the low-energy scale Mz (the Z-boson mass scale). 

The Standard Model (SM) suffers from the hierarchy problem associated with the need to keep electrowcak 
breaking much below the Planck scale. This problem is evaded if the theory is supersymmetric, with supcrsym- 
metry broken close to the electroweak scale. For this reason we consider the radiative corrections to neutrino 
masses and mixing in the context of the minimal supersymmetric extension of the Standard Model (MSSM) . 
We specify the low-energy boundary conditions of the renormalization group equations to be consistent with 
the three gauge coupling constants and the quark and lepton masses jl6| . We assume an effective SUSY scale 
of Ms = 500 GeV. We use the SM renormalization group equations below Ms and the MSSM renormalization 
group equations above Ms- The only boundary condition set at the family-symmetry breaking scale Mp is 
exact TBM mixing for the leptons []. The neutrino masses are set at the low-energy boundary relative to the 
lightest neutrino mass state (ra Vl with a normal hierarchy and m„ 3 with an inverted hierarchy). We keep 
|Am 2 2 | the solar mass difference and [Am^j the atmospheric mass difference. 

Figure |l| shows two contour plots. For the normal hierarchy the plot shows m Vl versus tan/3 and for the 
inverted hierarchy shows m„ 3 versus tan f3. The contours specify Log 10 Mp where Mp is the highest-energy 
family-symmetry breaking scale at which we can set TBM mixing and have the low-energy mixing angles 
consistent to within 4<j of the low-energy observations. The solar mixing angle #12 is the most sensitive to 
radiative corrections. Exact TBM mixing gives tan 2 £q 2 = 0.5, and our 4cr requirement at low energy translates 
to tan 2 6»i2 = 0.47 ±0.2 §. 

The difference between the two graphs can mostly be understood by the slight bias of the observational 
data (tan 2 612 < 0.5) and the opposite directions in which the tan 2 612 runs between the normal and inverted 
hierarchies. Starting with perfect TBM mixing at the scale Mp a normal hierarchy has tan 2 6*12 become 

1 We ignore the small departures from TBM at the high scale which may arise from diagonalising the chargcd-lcpton mass 
matrix [lj. 
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larger as the the renormalization scale becomes smaller. Once one falls below Ms, then tan 2 6*12 begins to 
get smaller as the renormalization scale goes down the during the final leg. The inverted hierarchy has the 
opposite behavior. Because there is a longer region of supersymmetric running, there is more parameter 
space of m V3 versus tan (3 compatible with Mp > 10 16 GcV. The slight bulge visible in the upper-right of the 
inverted-hierarchy contour plot with Mp « 10 4 GeV is due to the opposite directions of the supersymmetric 
running and the standard model running on tan 2 #12 in the region where Mp m Ms- 

The contours in Figure [l]hold implications for QD TBM family-symmetry models. For m vi > 0.1 eV, the 
neutrino spectrum is referred to as quasi-degenerate (QD) jl8j []. If cosmological observations are considered, 



they constrain the sum of the neutrinos JV ^ 0-42 eV at the 95% confidence level 19 . This implies 
m Vl < 0.14 eV which excludes the right half of Figure [j]. The remaining allowed narrow strip is consistent 
with the non-observation of ncutrinoless double beta decay (3(3q v which places a limit of m ee < 0.34 eV. 
Uncertainties in nuclear matrix clement weaken this bound by about a factor of 3. If we believe that the 
family-symmetry scale is greater than Mp > 10 10 GeV, and hypothesize a model which leads to a QD 
neutrino spectrum with normal hierarchy, then tan/3 < 6 (or tan/3 < 8 for a model with inverted hierarchy). 
In contrast, if a normal hierarchy model has tan/? > 6 (or tan/3 > 8 for an inverted hierarchy model), then 
the lightest neutrino need be less than 0.1 eV and therefore hierarchical. 



3 A discrete non-Abelian family symmetry model of QD neutrinos 
with TBM mixing 

As stressed in J20J an underlying 50(3) family symmetry readily leads to a near degenerate neutrino mass 
spectrum. In their model the chiral superfields, L l (where i is the 50(3) family index), contain the lepton 
doublets and transform as triplets under the 50(3) group. The chiral superfields containing the conjugates 
of the right-handed electron, muon and tau, respectively e c , [i c and r c , are 5*0(3) singlets. The effective 
Majorana neutrino mass is constrained by the symmetry and comes from the superpotcntial 

W eff =y (VL l )H u H u /M (1) 

where H u is the supermultiplet containing the Higgs field whose vacuum expectation value (VEV), (H u ) = v, 
is responsible for up quark masses in the MSSM and M is the messenger scale associated with the mechanism 
generating this dimension 5 term (in the Type II sec-saw it is the mass of the exchanged isotriplct Higgs field) . 

The important point to be taken from eq.(Q) is that the family symmetry forces the three light neutrinos 
to be degenerate. Small departures from degeneracy result when the 5*0(3) family symmetry is broken. In 
what follows we will show how this can naturally lead to a mass mixing matrix which gives near TBM mixing. 
This is done through the breaking of the family symmetry by the non-vanishing vacuum expectation values 
(VEVs) of familon fields, denoted as (f> l A , where the A = 3, 23, 123 labels three distinct fields and serves as a 
reminder of their VEV directions which are given by 





(2) 



where a, b and c are complex parameters. Table [l] lists the full set of supermultiplets and their symmetry 
properties under the 50(3) symmetry extended by a further set of symmetries G = Z^r x Z2 x U T (1) which 
limit the terms that can appear in the superpotential. Z3R is a discrete R— symmetry which ensures the 
familon fields arc moduli and cannot appear in the superpotcntial except coupled to "matter" fields carrying 
non-zero R— charge. The U T (1) symmetry is introduced to distinguish the third family of leptons from the 
first two. In practice it also explains why the mixing in the chargcd-lcpton sector is different from that in the 
neutrino sector which leads to near tri-bi-maximal mixing. 

The special structure of the VEVs in cq(||) is what will generate TBM mixing and is clearly the most 
important aspect of the model. This can happen naturally if the underlying family symmetry is not 50(3) 
but a discrete non-Abelian subgroup. We will discuss below the nature of this symmetry and the vacuum 



2 We define QD as m vi > 0.1 eV because above this value the /3/3rjK constraints for differing hierarchies and phases converge 
to a common region, as shown in figure El 
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Tabic 1: Assignment of the fields under the 50(3) family symmetry. 



alignment leading to eq(||) (the X field of Table |j] is introduced to facilitate this vacuum alignment), but first 
we show that it does generate approximate TBM mixing. 

The leading terms in the superpotential responsible for neutrino masses that are invariant under the family 
symmetries are given by 

W v = y {VV)H u H u + y Q (^\ 23 V) 2 H u H u + y @ {cj> 23 U) 2 H U H U . (3) 

where we have suppressed the messenger scale. Note that due to the Z 2 factor there are no cross terms 
involving </>230i23 |2l|, and due to the U T (1) factor there is no term involving ^3. As in eq(Jl|) , the QD 
mass scale is set by the first term of eq(||). For near degeneracy, the other terms must be relatively small 
(y©c 2 , y&b 2 -C yo, still suppressing the messenger scale). 
The chargcd-lcpton masses come from the superpotential 

W e = A e (£% 3 )e c #d + A^L'fe)^^ + K(L l cf> 3 )T c H d . (4) 

The m^/rrir ratio is given by \ij,{4>\^ l^d{4>\)- Using this the mixing between the second and third families 
of charged leptons is small of 0(m^/m T ). Similarly one may see that the mixing between the first and 
second families is of 0(m e /m ll ) and that between the first and third families is of 0(m e /m T ), both very 
small. Ignoring the small corrections from the charged-lcpton sector, the light neutrino mass eigenstates are 
proportional to the combinations 4>\ 23 L' l H u and <j) 23 L l H u 0. From cq(||) we see that these are given by 

v @ = — - v T ) (5) 
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where v e ^, T are the components of L e '^' T respectively (selected by the VEV of H u ). Ignoring the small 
chargcd-lepton mixings discussed above, v e .^.T can be identified with the current eigenstates. If b and c are 
real and positive, and to q = j/qc 2 w 2 < m@ = y@b 2 v 2 , one can see from eq(^) and eq(^) that we obtain the 
normal hierarchy, in which v@ may be identified with the atmospheric neutrino with bi-maximal mixing while 
Uq may be identified with the solar neutrino with tri-maximal mixing. The normal hierarchy persists for a 
range of complex b and c values in the neighbourhood of the real solution. An inverted hierarchy is possible 
and viable if 6, c are approximately imaginary and real, respectively. 

Although here we are working at the effective Lagrangian level, we already noted that (L l L l )HH naturally 
arises from the 50(3) invariant Type II see-saw mechanism. The other two neutrino mass terms can arise 
from Type I see-saw through exchange of appropriate heavy right-handed Majorana neutrinos, in a manner 
similarly to that discussed for a SU(3) based model in p3. Being of different origin it can readily happen 
that the common mass, toq = yov 2 is much larger than m@ and m,Q. 



3 In finding the mass eigenstates with a complex Majorana mass matrix, one needs to be careful to diagonalize M V M^ and not 
just Af„. Because M v is symmetric, it can also be diagonalizcd by an ortho gon al transformation OM u O T . In general O ^ U v 
and the square of the eigenvalues of M v are not the same as those of M v Ml & 
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Table 2: Action of the group factors Z 3 and Z 2 on the triplet representation 0\ 



4 Discrete non-Abelian symmetry and vacuum alignment 

We turn now to a discussion as to how the pattern of VEVs displayed in eq(|^) is dynamically generated. 
This can be achieved relatively simply if the underlying family symmetry is a discrete non-Abelian subgroup 
of SO(3) (and SU(3)). A very simple example is given by A4 = A(12), belonging to the A(3n 2 ) family of 
groups p5| . The A(12) invariant terms in the potential are those invariant under the group elements of the 
semi-direct product Z3 x Zi (which generate the group A (12)). The action of these group elements on a triplet 
representation i=1 ' 2 ' 3 is shown in Table ||. 

Since A(12) is a subgroup of 5*0(3), all 50(3) invariants are allowed by the discrete subgroup. Thus the 
terms of eq([|) and eq(Q) are allowed. The discrete subgroup allows additional terms, but these are all higher 
dimensional and consequently small provided the VEVs of eq(^) are small relative to the relevant messenger 
mass. Thus the lepton mass and mixing structure discussed is a consequence of the non-Abelian discrete group 
even though the 15*0(3) structure used above to motivate it is only approximate. 

Turning now to the question of vacuum alignment, consider the leading terms in the potential for the triplet 
familon fields. Because of the R— symmetry, in the absence of the X— field, there are no F— terms involving 
just the familon fields coming from the superpotcntial. The leading D— terms consistent with symmetries of 
Table are 



V(0) = am 2 ^2\c/) l \ 2 + /3m 2 



E 



7m 2 ^ 1 0' I + Sm 2 



(6) 



Here the quadratic term is driven by supersymmetry breaking and m is the gravitino mass. The coefficient 
includes radiative corrections which can drive it negative at some scale A, triggering a VEV for 0. The 
remaining terms can arise through radiative corrections and also only arise if supersymmetry is broken - hence 
the factor of m 2 on every term. The second term is generated at one-loop order if the superpotential contains 
a term of the form ^Yj^i^X 1 where \ % an d V are (Z^r = 1) massive chiral superfields which we take for 
presentational simplicity to have mass M. These two terms are invariant under the larger group 5C/(3) and, if 
a is negative, generate a VEV of the form (0) = (r, s, t) where r 2 + s 2 + 1 2 = x 2 , with x 2 a constant of 0(A 2 ). 
The third term, consistent with the non-Abelian family group, breaks SU(3) and 50(3). It will be generated 
if the underlying theory contains a superpotcntial term of the form ((0 1 ) 2 + w 2 (0 2 ) 2 + w(0 3 ) 2 )Z, where ui is 
the cube root of unity (ui 3 = 1) and Z is in a singlet representation of A(12) (one of the three irreducible 
singlet representations and distinct from the representation of Y) ^|. This coupling is invariant under the 
discrete group but not under SU{2>) or 50(3). The resulting third term of eq(||) splits the vacuum degeneracy. 
For negative a, the minimum for 7 positive has |(0 l }| = x(l, 1, 1)/V3 while for 7 negative l(0 l )| = #(0,0,1). 
Finally the fourth term also results from a one-loop radiative correction due to the £Y ^V 0*x* interaction. It 
is 5*0(3) but not SU(3) invariant and constrains the phases of the familon fields. For 5 negative and 7 positive 
the minimum has (0*) = 35(1, 1, 1)/V3 where x can be complex. This provides a mechanism to generate the 
vacuum alignment of 03 and 0123 as each will have a potential of the form in cq(|^) - as we are considering more 
than one familon, we label the coefficients with the familon's subscript to identify which term they correspond 
to. The structure of eq(|^) results if 73 is positive and 7123, ^123 are negative. 

Finally what about 023? Its VEV of the form in eq(Q) readily results once one includes the effect of the X 
field of Table |l| because the symmetries allow a term in the superpotential proportional to A(0230i23). This 
leads to a positive semi-definite term in the potential proportional to |0230i23| • There remains the need to 



align 023 and 03- This is readily done if radiative corrections generate a term m 2 



t 2 

6J, 3 03 



with a negative 



4 One may readily check that it is easy to assign charges under the group G to the new fields, xN Y an d Z to allow these 
couplings. 
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Figure 2: Neutrinoless double-beta decay mpp plots, from |l6|. From left to right, m m i n is the absolute 
value of the lightest neutrino mass, M is the sum of the light neutrino masses, and (mp) is the average mass 
determined from low energy beta decays. The shaded areas has width due to the unknown Majorana phases 
and the areas enclosed by solid lines take into account the errors of oscillation data. The two sets of solid lines 
correspond to the normal and inverted hierarchies. 

coefficient, thus a VEV will develop for ^23 in the direction given by eq(H|) []. One may readily check that the 
higher dimension terms allowed by the symmetry which involve the X field always involve an odd factor of 
023^123 and do not disturb the vacuum alignment mechanism discussed above. 

5 Neutrinoless double-beta decay 

The implication for neutrinoless double-beta decay in this model is unambiguous because the relative phases 
of the familon fields are determined. The amplitude for neutrinoless double-beta decay is proportional to the 
magnitude of ^2 m v% = "1/3/3 and this is what is measured. For TBM mixing U eT vanishes. The relative phase 
between the remaining two terms is given by Arg[m + e 2lpi2a itiq] — Arg[m ] where P123 = ^4r.9[yo</>i230i23/yo]- 
As m Q < to the relative phase remains small. This corresponds to the upper branches of Figure I in the QD 
region. Complex phases in the VEVs induce other CP violations through the charged-lepton sector that do 
not significantly affect m/3/3. 

6 Conclusion 

Attempts to explain the structure of fermion masses and mixings often rely on structure at a high (Grand 
Unified?) scale, Mp, to generate the observed pattern. One possibility, consistent with neutrino oscillation, 
is that neutrinos are nearly degenerate. However, due to enhanced radiative corrections in this case, the 
observation of near TBM mixing is difficult to reconcile with such a high scale mechanism. To keep the 
deviations from TBM mixing within experimental limits it is necessary to limit the scale at which TBM 
mixing is generated. We have determined this scale for the MSSM and found significant bounds on Mp. For 
example, for degenerate (m Vi > 0.1 eV), normal hierarchy neutrinos and tan/3 > 6 (or tan/3 > 8, for inverted 
hierarchy neutrinos), Mp < 10 10 GeV is required. To get close to the Grand Unified scale with QD neutrinos 
it is necessary to have very small tan f3. 

5 Noto that in eq(^) wo have used the freedom to define the directions such that {4>\' 2 ) = 0, (^23) = 0- 
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Turning to the origin of the structure, we have constructed a model based on a discrete non-Abelian family 
symmetry which gives a QD neutrino spectrum and near TBM mixing. This relics on a natural mechanism for 
vacuum alignment of the familons which break the family symmetry. The mechanism predicts that neutrinolcss 
double-beta decay should be maximal. Although we only constructed the low-energy effective theory, it fits 
very well with a see-saw mechanism in which the degenerate mass comes from a Type II see-saw while the 
small departures from degeneracy are driven by a Type I see-saw. 
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